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It is shown that a 25(20)% difference between the decay constants /d 3 (/b s ) and 

/d(/b) occurs due to large differences in the pole masses of the s and d(u) quarks. 

The values rjo = foj fo ~ 1.23(15), recently observed in the CLEO experiment, 

' and rjB = fsJfB ~ 1-20, obtained in unquenched lattice QCD, can be reached 

G\' 

■ only if in the relativistic Hamiltonian the running mass m s at low scale is m s (~ 0.5 

(N 

GeV)= 170 - 200 MeV. Our results follow from the analytical expression for the 
pseudoscalar decay constant fp based on the path-integral representation of the 
• meson Green's function. 

a: 

Oh! 

Relativistic potential models (RPM) have been successful in their description of light- 



light and heavy-light (HL) meson spectra [1, 2]. Still there exists a fundamental problem, 



^ . which remains partly unsolved up to now. It concerns the choice of the quark masses in the 
kinetic term of the relativistic Hamiltonian, which in different RPMs vary in a wide range. 
For example, HL mesons were studied with the use of the Dirac equation, taking for the light 
quark mass m n (n — u, d) — 7 MeV in [3] and 72 MeV in [4], and in the Salpeter equation for 
the strange quark mass the values m s = 419 MeV in [5] and m s = 180 MeV in [6] have been 
used. However, in contrast to constituent quark models, where the constituent mass can 
be considered as a fitting parameter, a fundamental relativistic Hamiltonian has to contain 
only conventional quark masses-the pole masses. These masses are now well established 
for heavy and light quarks [7]. They have been used in the QCD string approach giving a 
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good description of meson spectra [6, 8, 9]. However, the strange quark mass m s is still not 
determined at low scale. At present, owing to the QCD sum rules calculations [10] and lattice 
QCD [11], m s is well established at a rather large scale: m s (/x = 2 GeV) = 90 ± 10 MeV, 
while in the Hamiltonian approach the mass m s enters at a lower scale, which is evidently 
smaller that the scale /i c 1.2 GeV for the c quark. Therefore it is very important to find 
physical quantities which are very sensitive to m s at low scale: fj, s < 1 GeV. In this letter 
we show that such information can be extracted from the analysis of the decay constants of 
HL mesons, namely, from the ratios JdJ fn and f B J f B . 

Recently, direct measurements of the leptonic decay constants in the processes D{D S ) — > 
fiu^ [12, 13], and B — > Tfi T [14, 15] have been reported. In Refs. [12, 13] the CLEO collab- 
oration gives f Ds = 274(20) MeV and f D = 222.6(20) MeV with r] D = f D Jf D = 1.23(15), 
having reached an accuracy much better than in previous experiments [16]. This central 
value for i] D appears to be larger than in many theoretical predictions which typically lie in 
the range 1.0 — 1.15 [17, 18, 19, 20]. Therefore, one can expect that precise measurements 
of rjD and tjb in the future can become a very important criterium to distinguish different 
theoretical models and check their accuracy. In particular, relatively large values 

VD = 1.25(3), rj B = 1.20(3), (1) 

have been obtained recently in lattice (unquenched) calculations [21, 22] and also in our 
paper [6]. In this letter we show that 

1. The running mass m s {^i) at a low scale, /ii ~ 0.5 GeV, can be extracted from the 
values i] D and i] B , if they are known with high accuracy, < 5%. 

2. The values rjo and as given in Eq. (1), can be obtained only if the running mass 
m s (/xi) lies inside the range 170 — 200 MeV. In particular, in the chiral limit, rrid = m u = 0, 
as well as for m d = 8 MeV, and for m 3 (fXi) = 180 MeV, the ratios r] D and t]b calculated here 
are 

^ = 1.25, 77b = 1-19. (2) 

3. The value m s (0.5 GeV) satisfies the relation 

m s (0.5 GeV) 1 nn 

\ - — i » 1.97. 3 
m s (2 GeV) V ; 

In our analysis we use the analytical expression for the leptonic decay constant in the 
pseudoscalar (P) channel, derived in Ref. [6] with the use of the path-integral representation 
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for the correlator G P of the currents jp(x): G P (x) = (jp(a;)ip(0)) vac : 

J P = / G P (x)d 3 x = 2iV c V (yp) " lv? " (0)|2 e- M " T , (4) 

where M n and tp n {r) are the eigenvalues and eigenfunctions of the relativistic string Hamil- 
tonian [23, 24, 25], while oJ qn {ojQ n ) is the average kinetic energy of a quark q(Q) for a given 
nS 1 state: 

Uqn = (^m 2 q +p 2 ) nS , UJ Qn = (^m 2 Q +p 2 ) nS . (5) 

In Eq. (5), m q {mQ) is the pole mass of the lighter (heavier) quark in a heavy-light meson. 
The matrix element (lp)„ refers to the P channel (with exception of the n and K mesons 
where additional chiral terms occur) and was calculated in Ref. [6], 

(Y P ) n = m q m Q + UJqnUQn - (p 2 )nS- (6) 

On the other hand, for the integral Jp (4) one can use the conventional spectral decompo- 
sition: 

Jr= [ GA-)d 3 x = 2^-(/ P ") 2 e- M " T (7) 

J n n 

Then from Eqs. (4) and (7) one obtains that 

2JV c (y P ) B |y n (0)| 2 

All necessary factors in Eq. (8) for the ground state (n = 1) and the first radial excitation 
(n = 2) are calculated in Ref. [6] but here we consider only ground states and omit the index 
n everywhere. Our calculations are performed with the static potential Vo(r) — or — | aB ^ 
[9], and the hyperfine and self-energy contributions are considered as a perturbation. It is 
important that our input parameters contain only fundamental values: the string tension a, 
the QCD constant A(nf = 3) in a B (0) an d the conventional pole quark masses. For the s- 
quark mass m s {^\) one can expect that the scale /xi is close to the characteristic momentum 
~ V (p 2 ) ~ 0-5 — 0.6 GeV. This scale also corresponds to the r.m.s. radii Rm(1S) of the 
meson we consider. For HL mesons 

R D « R Ds = 0.55(1) fm, R B ss R Bs = 0.50(1) fm, (9) 

so that iii ~ R^} 0.4 — 0.5 GeV. We show here that this mass m s (/ii) is strongly correlated 
with the values of rj D and r] B . For other quarks we take m c = 1.40 GeV and m b = 4.78 GeV 
[8]. 
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It is of interest to notice that for HL mesons the ratios 



to = Zo s = = 0.347(3) (10) 



UJqUJc 



are equal for the D and D s mesons with an accuracy better than 1%, and also that these 
fractions for B and B s mesons coincide with < 2% accuracy (vf(0) = i?f(0)/47r): 

& = U = = 0.146(2). (11) 

UJqUJb 

It is important that the equalities £d = £d s and = £b s practically do not depend on the 
details of the interaction in HL mesons. Therefore, in the ratios i]d{vb) the factors given in 
Eq. (10), £d(£b) cancel and one obtains 

2 _ / m s m c(b) t u s uj c{b) - (p 2 )d s{Bs ) \ M d(b) 

Vd(b) ~\(y p } d{b) + (YpW) J m Ds{Bs) - {U) 

In Eq. (12) the second term is close to 1.05, while the first term, proportional to m s , is 
not small, giving about 30-60% for different m s (below we take m s from the range 140 ± 60 
MeV/c 2 ). With an accuracy of < 2% 

rf D = 2.708 x m s (GeV) + 1.07(1), if m d = m u = 0, 
rf D+ = 2.648 x m s (GeV) + 1.05(1), if m d = 8 MeV, (13) 

i.e., in the chiral limit 

t)d — 1-14 (m s — 85 MeV), 1.25 (m s — 180 MeV), 1.27 (m s — 200 MeV), (14) 

and for = 8 MeV, t\d = 1.13, 1.24, and 1.26, respectively, for the same values of m s , so 
decreasing only by ~ 1%. 
For the B and B s mesons 

rj 2 B = 1.90 x m s + 1.07(1) (m d = m = 0); 

rf B0 = 1.871 x m s + 1.07(1) (m d = 8 MeV), (15) 
which practically coincide, and in the chiral limit (m d = m u = 0) 

r] B = 1.11 (m, = 85 MeV), 1.19 (m s = 180 MeV) 1.21 (m s = 200 MeV). (16) 
These values of r] B appear to be only 3 — 5% smaller than r\ D . 
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Thus for m s = 180 MeV and m d = 8 MeV we have obtained 



VD+ = 1.25(2), r] B = 1.19(1), 



(17) 



in good agreement with the CLEO data: i] D (exp) = 1.23(15) [13]. 

To check our choice of m s = 180 MeV, we estimate the ratio m s (0.5 GeV)/m s (2 GeV) 
using the conventional perturbative (one-loop) formula for the running mass [26] 

-d„ 



m{fi 2 ) = m ( ^ In ^~ 



lnlnfe 
In 



2~ A 2 

Here m is an integration constant and the other constants are 



(18) 



19 \ 2 
51 - Y n f J > A) = 11 - d r, 



4 

Jo' 



(19) 



To calculate m s (2 GeV) we take rif = 4, A(nf = 4) = 250 MeV (/3 = y, d m - 
di = 0.3548). We take n/ = 3, A{n f = 3) = 280 MeV to determine m s (l GeV) and m s (0.5 
GeV) (for n f = 3, /3 = 9, rf m = |, d x = 0.3512). Then, from Eqs. (18,19) m s (2 GeV) 
= 0.618 77io, m s (l GeV) = 0.7825 m , and m s (0.5 GeV) = 1.217 m and therefore 



m s (l GeV) 



1.27, 



m s (0.5 GeV) 



1.97 . 



(20) 



m s (2GeV) ' m s (2 GeV) 

It means that m s (0.5 GeV) = 180 MeV, which we used in our calculations, corresponds to 
m s (2 GeV) = 91 MeV which coincides with the conventional value of m s (2 GeV) = 90±15 
MeV [7]. Thus our estimate of m s (0.5 GeV) = 180 MeV supports our choice of this value 
in the relativistic string Hamiltonian, which provides a good description of the HL meson 
spectra and decay constants, and gives rise to the relatively large values of r\ D and r\ B in 
Eq. (1). 
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